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1.0 INTRODUCTION 


The ocean tides, with 'their principal semidiurnal and diurnal 
periods and their varying geometric patterns, act as sources for de- 
forming the earth's crust and mantle. Understanding these deformations 
has importance for geodesy, geodynamics and oceanography. Moreover, 
information about the interior of the earth may be obtained through the 
magnitude and frequency dependence of the response of the solid earth 
to the tidal forces. The tide of the solid earth is composed of 
two basic responses; (1) a body tide due to the yielding of the earth 
to direct forces of the sun and moon and (2) a "load" tide produced by 
surface forces from the time varying ocean tides. It is difficult to 
distinguish between these two responses because their time dependence 
is similar, being driven by the same ultimate force. However, the nat- 
ure of the driving force of the "body tide" is well understood, while 
the knowledge of the deep ocean tides through global numerical modeling 
is a recent advancement (Perkins and Accad (1969), Hendershott (1972), 
Zahel (1977), Estes (1977), Schwiderski (1978), Parke (1978)). The body 
tide varies in a relatively smooth nature over the earth's surface, de- 
pending principally on averaged overall elastic properties while the 
load tide is complicated by discontinuities of the surface load at 
coastal boundaries and by local ocean tide circulation (e.g. amphi- 
dromes and anti -amphi dromes ) . Moreover, the displacement of the load 
tide is appreciable only in the crust and upper mantle, while the body 
tide has relatively large amplitude through most of the earth's 
interior. The load tide response then depends more on local crustal 
properties so that variations in near surface earth structure will be 
more reflected in the load tide, 

The response of the solid earth to ocean loading may be evaluated 
by convolution of a model for the ocean tide over the global ocean 
with appropriate Green's functions which are derived from models of 
the earth's interior. In the present study, ocean tide models are 
used together with the Green's functions calculated by Farrell (1972) 


In terms of a layered spherically symmetric Gutenberg-Bull en earth 
model to calculate global values of horizontal and vertical crustal 
displacement, gravity perturbation and strain at the earth's surface. 

Recent global ocean tide models differ in detail and the accuracy 
with which the tides may be predicted from these models is uncertain. 
However, a measure of a model's overall accuracy may be provided by 
its agreement with recent astronomical calculations for the lunar 
acceleration and the rate of energy loss, which are simply related to 
the amplitude and phase of the (2,2) tesseral harmonic of the tidal 
elevation. Goad and Douglas (1977) have analyzed perturbations in 
satellite orbits which are also proportional to the low order harmonics 
and obtained values which show very close agreement with the astronomical 
values and the values calculated from the Mg models of Schwiderski 
(1978) and Estes (1977). Although this provides some confidence in 
these two numerical Mg models, it must be pointed out that close 
agreement with the (2,2) harmonic in a tidal model does not necessarily 
imply a correct model of the tides in specific regions. Tide solutions 
computed from Laplace's Tidal Equations fall into two groupsj those 
constrained to agree with coastal observations and those which employ 
no data or constraints. The Estes (1977) tide models selected for 
the present study are of the second category and provide solutions 
for the Mg, Sg, Ng, Kg, K-j , O-j and constituents integrated at 
3° special resolution. The effects of ocean loading and self-gravitation 
have been included ia deriving the Mg model. These models are 
reasonable candidates for the global calculations of earth response 
for inland regions and open ocean areas, which are sensitive to the 
large scale effects of mid ocean tides. However, for coastal areas 
where the response is strongly influenced by local water tides, a 
model which provides a finer special resolution and incorporates coastal 
data constraints, such as a regional empirical model or the global 
Schwiderki (1978) model will provide greater accuracy to the convolution 
computation. 
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The results of the global calculations for crustal displacements, 
gravity anomalies and strains caused by ocean loading are presented 
in the form of corange and cotidal charts. The phases described by 
the cotidal contours are relative to the Greenwich meridian, and are 
expressed in hours instead of angular measurement. Hour values are 
obtained by dividing the phase expressed in degrees by 15 degrees 
per solar hour. In addition, a software package has been developed 
which will evaluate the vertical displacement due to loading by the 
principal tidal constituents and the solid earth tide as a function 
of geographic position and time as specified by user imput, 


2.0 EARTH RESPONSE TO OCEAN LOADING 


The calculation of the earth deformation due to surface mass loads 
closely follows the Green's function approac^/of Longman (1964) and 
Farrall (1972). Farrell hds integrated the equations of motion for a 
self-gravitating elastic spherical earth using a Gutenberg-Bull en A 
earth model and produced load Love numbers h|j , 1' , and k' to 
high order n . The elastic earth response then reduces to a convolution 
of the ocean tide with the areed's function over the global ocean, ' 

s(+.x;t) • |(r^ pGfWjdd' (1) 


where 


cosf ® + aos^cos^'co8{\""^ ) 

and 5 ((j),x;t) denotes the ocean tide and R is the mean radius of 
the earth. Here G^(y) represents the point load Green's function 
transformed to an earth-fixed coordinate system and R represents the 
appropriate response. The total tide is approximated as a sum of con- 
stituent tides 

e((^,X;t) I 5 . (*,X;t) ( 2 ) 

where 

e^(4»^5t) = i4^{f,x) c£7s[o.(t-t^)-\|i^.(<}),x)+c^] (3) 

and constituents i = {Mg, Sg, Ng, Kg, Kp 0^ } are available 

from global numerical models. Here t^ is January 0, 1900 and the 
phase function ip is relative to the lunar passage at Greenwich of the > 

/ 
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Table I 

Constants for Constituent Tide Models 










the ficticious moon for the particular constituent. The frequencies 
ami constant phases are presented in Table I. The total 
response is then 


1 ' 

where 


( 4 ) 


i?,-(4.X;t) » oo8[a^.(t-tQ)+e|]+Q^ (t-t^l+e^] (5) 

Gf G, 


and 




( 6 ) 


In terms of amplitude and phase , the response to the i 
constituent is 


ith 


i?^(>,A;t) - A^{ij),X) (3oe[cJ^(t-tQ) - + e^] 


( 7 ) 


v^here 


A,(*.*) =/ Pf + 


Gf ’Gf 


■« 1 , 


A^{(j>»A) = tan 


-1 


( 8 ) 
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The appropriate Green's functions for the augmented potential and 
surface vertical and horizontal displacements at an angular distance 
Y (spherical earth) from a point load at the pole per unit of loading 
mass are 


*"(y) ' ^ I k;; P^Ic^’sy) 
% n=0 " " 


U"{y) « if I h' P (aos-i) 
e n*=0 

v'(y) “ if I ir 


e n=l 


n 3 y 


C9) 


where is the mass of the earth, R is the mean earth radius and 
g is the acceleration of gravity at the surface. Here primes on the 
Green's function denote that they are in a symmetric point load coordinate 
system. The Green's function for the differential gravity acceleration 
as given by Farrell (1972) is 


G'(y) = Jo (n+Zh^-(n+1)k') P„(ooay) . (10) 

Here 6"(y) represents the difference between g , the acceleration of 
gravity at the earth's surface, and the acceleration on the deformed 
surface after application of the ocean load. Following Farrell, we 
break the acceleration Green's function into direct, or Newtonian 
acceleration 

= if I n PA<^osy) (11) 

e n=0 " 


and elastic acceleration 


( 12 ) 


G'^t) . G'(r)-G'N(T) 

■ (?- Z C2hj;-(n+l)k'] p (oosf) . 
e n*0 


As pointed out by Pekeris (1978), equation (10) ap;jl1es to gravity 
measurements made at a point below the tidal loading sheet. For points 
above the sea, the increase In gravity caused by crossing the loading 
sheet from below 

-4ir G p 5 (13) 


must be added to equation {1$) which results In 

SMy)-#- I C2h'-(n+l)(1+k')] , (14) 

i^e n=0 " ^ 


The correction term of equation (13) may be written 

" g ^ (15) 

P ^ 


where p denotes the density of sea water and p the mean density 
of the earth. For gravity measurements on the coast this correction 
is important. However, for inland measurements where there is no 
ocean tide the correction term is zero. Pekeris (1978) also points 
out that in evaluating the Newtonian acceleration Green’s function of 
equation (11), Ferrell (1972) has omitted the delta function from the 
expression 



n P^(x) - 


1 


+ 6(l-x) 


and thus a term 


||.g| ■ (16) 

should be added to the gravity perturbation. Again, this term con- 
tributes only over ocean areas. Similarly, expressions may be obtained 
for the Green's functions of the non-zero elements of the strain tensor 
at the earth's surface. As polntedmtby Farrell ( 1972), these 
Green's functions are slowly convergent series and must be summed to 
large values of n , 


The summed values for the Green's functions as a function of 
angle have been taken from tables given by Farrell (1972), where 
values are available for U'', V^, 6' and , and. S" is the 

strain tensor. The other diagonal components of the strain tensor 
(off-diagonals are zero for the syumetric point load coordinate 
system, with load at the pole, for which the Green's functions are 
derived) are calculated from 



(17) 


where x(a) and o(a) are the Lame parameters at the top layer of 
the earth model. Here y is colatitude and X longitude in the 
symmetric point load coordinates. To evaluate the convolutions for 
the earth response of equation (1), it must be realized that the 
primed Green's functions described above are with respect to the 
symmetric point load coordinate system, and to resolve components of 
vector and tensor quantities appropriate transformations must be 
applied. Let (<|>,x) be the latitude and longitude of the point of 
evaluation for the convolution and be the latitude and 
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longitude of the water column being considered as the load. Moreover, 
let r , 0 , X denote the unit vectors In the direction of Increasing 
r, 0, X at the point of evaluation (where e Is the co-latitude), 
Then the Green's functions for the horizontal displacement In the e 
and X directions, respectively, are 


Ve(T) “ V‘'(if)coao 
Vj^(y) * V'(y)Bina 


( 18 ) 


where 


-- _ stnit'-sinicoBr 



eina « 

6%ny 


(19) 


and as before 


cosy = ein^sin^' + co8^co3ycos(X-X') 


(20) 

einy ~ /l-cos^y 

Similar expressions for the radial displacement and strain tensor components 
are 


U(y) • U'(r) 

S ■ S' ■ 
rr rP' 

(21) 

See ■ 

Sxe ■ 

*xx * • 


The computation for the response convolutions of equation (1) is 
performed at each 3®x3® equal -angular grid point over the surface of 
the globe using the Estes (1977) 3® resolution models. Cotidal and 
eorange charts for the principal constituents are presented in Figures 
1-7. Numerically, the quadrature is evaluated as the sum of near zone 
and far zone contributions. The far zone is defined as the collective 
region of 3®x3“ surface areas whose centers are greater than 3° from 
the point of evaluation. Since the tide models define values only at 
equal-angular 3°x3® grid points, the tide amplitude and phase are 
considered constant within each 3°x3® surface area element. The far 
zone contributions to equation (1) are then obtained by summing the 
individual area elements with the Green's function evaluated at the 
central angle between the center of the element and the point of 
evaluation. In the near zone, the Green's functions change so rapidly 
that the elements of surface area are evaluated on a graduated scale 
from yjg- “ x to x ® as the integration proceeds awa/ 
from the point of evaluation. Again, the tide amplitude and phase 
are considered constant within each 3®x3° surface area. 

The computation for the earth response is greatly influenced y 
the near zone contribution in most ocean and coastal areas. In e 
open ocean areas the Estes (1977) tide models and the assumptic o-' 
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constant tidal amplitude and phase within a 3®x3® surface apta 
are adequate. However* near coastal zones the calculations suffer 
from the fact that the global numerical tide models used are theoretical 
solutions of the Laplace Tidal Equations for the open ocean areas of 
depths greater than approximately 500 meters and the solutions contain 
no data. Moreover, the 3® approximation to coastline geometry Is 
course, Accurate near zone values for coastal regions should be ob- 
tained from more detailed regional tide models or from finer special 
resolution global models, such as the M 2 model developed by 
Schwiderski (1978), that have incorporated coastal tide observations 
into the solution to more precisely represent the real near shore 
variation. Because of these considerations, the calculations pre- 
sented in this study for regions within 1® to 3® of coastlines should 
be viewed qualitatively in the transition from open ocean to inland 
areas. Note also that the tide models used in the calculations pre- 
dict only to 81® north latitude, and the tide is assumed to be zero 
above this limit. The responses evaluated in the north polar regions 
then are only the far zone contributions, Near zone contributions C 
should be computed from a special pole tide nwdel. 

2 . 1 Crustal Displacements 

The vertical response due to ocean loading by constituent “i" 
as given by equation (7) is 

U^(l.,X;t) » Ay((j>i.A) ae»sC(j^(t-tQ) + 

where the amplitudes and phases are calculated by convolution with the 
Green's function from equation (9) 

6f(Y) H U^y) 

Corange and cotidal maps for constituents M 2 , Sg, Ng, Kg, Kp 0i and 
P.j are presented in Figures 8 through 14, As indicated previously, the 
phases are with respect to lunar passage of Greenwich. Horizontal 
responses 
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Vx ihMt) 

calculated from the 6reen*s functions of equation (18) were evaluated 
only for the Mg tide and are presented in Figures 15 and 16. 

The vertical response exhibits a spacial behavior over the oceans 
which resembles the structure of the constituent tidal heights. The 
greatest amplitudes are on the order of 5 cm. in the region of the mid- 
Pacific anti-amphidrome predicted by the Mg model, while over con- 
tinental areas the calculated load response is generally less than 1 cm, 
and decreases with distance from the coasts. Horizontal displacements 
are an order of magnitude smaller than the vertical displacements and 
exhibit a more complex structure, 

2.2 Gravity Perturbations 

The measurement of gravitational acceleration is the resultant 
of components 

g « g^ + 9g + gg + g^ (22) 

where represents the acceleration on a rigid earth, gg represents 
the contribution from a symnetric, ocaanless, elastic earth, gg 
represents the contribution from oce!'.n tides, and g^ is the response 
of regional anomalous geologic structure. The components g^ and gg 
are the larger terms and have essentially the same phase, while the 
terms g^ and g^ may differ considerably from other terms with 
respect to phase. 

Modern gravimeters are caoable of detecting the g^ and g^ con- 
tributions with limited precision, and it is now clearly established 
that earth tidal gravity parameters for both coastal and inland measure- 
ments are sensitive to the regional and global ocean tides. A 
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comparison of accurately measured patterns of tidal gravity spacial 
variation within continental regions with model calculations will 
offer the ability to distinguish between tide models and establish 
their levels of precision. In fact, Kuo, et al. (1970) attempt to 
indirectly map the ocean tides by solving the inverse problem of the 
response of tidal gravity to ocean tide loading. The accuracy of 
data will be greatly improved with the use of a new type of super- 
conducting gravimeter (Warburton, r»t ;sl (1975)). Presently measure- 
ments using this instrument are available only at La Jolla and Pinon 
Flat, California. 

Calculations of tidal gravity at specified stations using 
available tide models and the comparison of results with observations 
have been performed by several authors. In particular, Robinson (1974) 
compared the observed relative gravimetric factors of M 2 and 0^ at 
several stations in the southeastern United States using five different 
published global Mo models. Wilson (1978) has extended the analysis 

ii. 

using additional ocean tide models and a model for the Gulf of Mexico. 
Bretreger and Mather (1978) have analysed tidal gravity measurements 
in Australia using a ten-parameter response model and global Mg tide 
models. However, only Kagan and Polyakov (1977) using an Mg model 
by Gordeyev, et al (1976) have performed calculations over the entire 
earth to present a global picture of variation. 

* ' 1 

The amplitude of the gravimetric factor 

6 = l- '^k + h 

where k and h are Love numbers, is the ratio of tidal gravity at 
a point on the earth to the theoretical amplitude on a rigid earth. 

Th?v local epoch of a measurement is the phase difference between the 
observed gravity and theoretical gravity. In the absence of ocean 
loading the value of the gravimetric factor would be approximatly 1.16. 
The correction as is defined by vectorial addition of the ocean 
load gravity vector to the theoretical gravity 
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/■'v 

c ] 


" y^Ocean load gravity 

/phase of ocean tidal gravity 
^ Local Epoch 

Oceanless Theoretical gravity 

where the g^ contribution is assumed negligible for gravity and 
the theoretical gravity has zero phase lag. Note that the phase of the 
ocean tidal gravity in this diagram is with respect to the phase of 
the theoretical gravity. The relation between the local phase of . 

theoretical gravity and the Greenwich phase is /r 

)) 

Local Phase - Greenwich Phase + mx 

where x is the east longitude of the measurement point and m is the 
order of the tide (2 for semidiurnal and 1 for diurnal tides). The ampli- 
tude of the ocean load gravity in ugal is obtained from the amplitude 
|a 6| multiplied by the amplitude of the theoretical gravity of the rigid 
earth. 


The ocean loading gravity perturbations for the M2 and 0^ tides 
are calculated from the Green's function of equations (11) and (12) 
for the Nev/tonian and elastic contributions and are presented in 
figures 17 through 20. Here the calculation for the Newtonian term 
neglects the delta function pointed out by P.ekeris (1978). For coastal 
and ocean areas the additional term given by equation (16) may be 
evaluated from the ocean tide maps given in Figures 1 through 7. 
Moreover, for measurements evaluated at coastal and ocean points 
above the tide sheet, the correction term given by equation (15) must 
be added to obtain the total gravity perturbations. The amplitude 
and. phases of the gravity perturbations, as with the horizontal 
crustal deviations, exhibits a spacial behavior which resembles the 
tide structure. This supports the suggestion that the complexity of 
the ocean tidal gravity is due to the complexity of the global tide 
rather than the response of the earth to the ocean tidal loading. 
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The responses are on the order of one ugal over continental regions 
and decrease with distance from the coasts. Values measured by 
Warburton, et. al. (1975) at Pifion flat for the 0^ and Mg tidal 
gravity components using the superconducting gravimeter are In 
reasonable agreement with the calculations presented In Figures 17 
through 20. However, their measurements on the coast at La Jolla 
(on the order of 100 km. distance from Pifion Flat) are substantially 
greater than the results shown in Figures 17 through 2p. This is 
most likely a consequence of the Inadequate precision In modeling the 
near zone ocean tide at the coasts as discussed In Section 2.0. 

2.3 Strains 

A comparison of precise tidal strain observations with model 
predictions could provide insight into the nature of the earth and 
its structure. However, obtaining accurate measurements of tidal 
strains is considerably more complicated than for tidal gravity. 

While g^ from equation (22) is considered negligible for tidal 
gravity, tidal strain measurements are very sensitive to local topo- 
graphic and geologic anomolous influences. A basic understanding of the 
ocean tidal loading strains and the local influences could have a great 
impact on the field of techtonic geophysics. In particular, Young and 
Ziirn (1979) claim to have provided weak evidence that earthquakes in the 
Swabian Jura are triggered by tidal shear stress. The concept that 
ocean load tide strains could provide a mechanism for triggering is 
feasible, as ocean loading can introduce appreciable horizontal shear 
strains when loading is not laterally uniform, such as near a coastline. 

The ocean leading tidal strain tensor have been evaluated for the Mg 
tide from convolution with the strain component Green's functions of 
equation (21). The symmetric tensor is written 

(23) 


S = 


e 0 

YY 

0 e 


0 

60 


•X0 '"XX 
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where each element varies with time as 

e. .(^»A;t) * A (♦jA) ooeCou (t-t )-A (4>>A)'^e|j ] . 

Note that the tensor has only three independent components at the free 
earth surface, since 



11 } < 




) 


where x(a) and o(a) are Lame parameters at earth radius, a . 
The horizontal surface strain tensor may be written in the form 


®ee ®ex 


' it®8eS.x> “ 

®xe ®AA 

\ J 


“ *<'ee*«xx> , 


^^®ee"®AA^ 

®ex 


®ex 

"^^®ee“®xx^, 

// 

(24) 

o(a) ^ r 

nftrv ' 

' i<^6-«xx> «ex 

+ 

®9A “^^®00"®XX^ 



as the sum of a pure areal strain and a pure shear strain. While it is 
conventional to represent the three independent components of the surface 
strain tensor by plotting the linear strain at a particular geographic 
position as a function of azimuth on a polar diagram, we present more 
detailed corange and cotidal plots of the elements e„„ , e*. , e,„ 

and in Figures 21 through 28. Due to plotting difficulties the 

tensor components are presented in separate figures over ocean areas 
and land areas. 
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The radial strain component shown in Figures 21 indicates 
strains of the order of 5x10”^® in the mid-oceans and IxlO"^® over 
continents. The transitions over, coastal areas appear generally smooth. 

This is in contrast to the e^g and components of Figures 22 

through 25. Here the mid-ocean strains are on the order of 5x10“® while 
strains over continental regions are two orders of magnitude smaller. 

The steep gradients occur at the shorelines, where abrupt changes in phase 
are also observed. Mote that the amplitudes of Ogg and are nearly 

equal while their phases are nearly 180® apart. This is consistent with e^^ , 
which is proportional to their sum, being two orders of magnitude smaller. 


e.«‘ strains of Figures 26 and 27 show a comparable magnitude to 

Aw 


and e^^^^ over continental areas, while being approximately an order 
of magnitude smaller over mid oceans. Figure 28 displays that the e^^^ 


strain generally shows a sharp increase in magnitude to the order of 
-8 

5x10" near coastlines before rapidly falling to smaller values over 


land areas-. Abrupt changes in phase also occur near shorelines. By 


equation (24), these values show that the ocean load strains are 
dominated by pure shear, while the areal strains are small; 



3.0 COMPUTER SOFTWARE FOR BODY AND LOAD TIDE RADIAL DISPLACEMENT 


A computer program has been developed to evaluate the radial dis- 
placements given in Section 2.1 due to ocean loading and the total 
radial displacement due to the solid earth tides (body tides) at a 
user specified position on the earth for a desired time interval. The 
'algorithm evaluates the contributions from the ocean tidal constituents 
and the body tide individually for a one day period at one hour incre- 
ments. Required input are the Modified Julian Day and the geodetic 
coordinates of the point of interest. 


The radial ocean loading displacements of Section 2.1 are represented 
in the software In the form of spherical harmonic expansions 


U^(i>,x;t) = Ay((^,x) cos[o^(t-tQ)+c^--A^( 4 »x )3 

* { I Ca„^coemX+b^jj^8inmx]Pj^jj^(ein4.)}cos(cr^(t-tQ)+e^.) 


n,m 


+ { I [Cj,^£JoemX+d^ij^8inmx]P (stn4»)}8in(a^(t-tQ)+e^) 

n,m 


(25) 


where the coefficients a^^, b^^^, c^^ and are obtained by a least 

squares fit to degree and order twelve. These expansions provide an 
accurate and compact method for evaluating the ocean loading tides. 

The body tide elevation is given by 


M^R 


UB(*,X-.t) = 4^ J h 


n=2 


R. 


n+1 


P^(CO80| 


MS^ 


(26) 


where M^j is the mass of the disturbing body (Moon or Sun), R^j is the geo- 
centric distance to the body and h^^ are Love numbers. The angle 
denotes the geocentric zenith angle of the moon (sun) from the point of eleva 
tion. The terms in the expansion fall off rapidly so only the first term 
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Is of major significance. The radial body tide, Ug , is then 
evaluated as 


Ug (^x;t) 






where r represents the unit radius vector at the point of Interest on the 
earth: 

r ® [cos^cosX, coe^sinX, ein^"] 


where geodetic latitude, ^ , and longitude, X , are program Inputs. 

represents the unit vector from the center of the earth In the 
direction of the disturbing body: 

Rjj " Cx ** » V*] 

where X", u", v" give the position in earth-fixed coordinates. 

In calculating R for the moon, a true longitude and the latitude 
(above the plane of the ecliptic) are derived from the Hill-Brown theory 
using the Modified Julian Day. Brown's tables express the coordinates of 
the moon as sums of periodic terms whose arguments are algebraic sums of the 
multiples of 1, i't F, D, r . See Tables II and III. 
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TABLE II 

Ecliptic Eletnents 


MJD ■ Modified Julian Day 
D » MJD-241 5020.0 
01 - D*l.E-4 

A « 296. 104608+13- 0649924465*D+0.0006889*(D1)^ 
I' « 358.475845+0, 9856002670*0-0, OOOOn2+(Dl)^ 
F » 11. 250889+13. 2293504490*0-0. 0002407*(Dlj^ 

D * 350. 737486+12. 1907491914*0-0, 0001076*(D1)^ 
r « 281. 220833+0. 470684*01+0. 339E-4*(D1)^ 
c « 23. 452294-0. 0035626*01-0. 123E-6*{D1)^ 

0g » 99.6904833+360.98564733*0-180.0 


Table III Development of Lunar Position 


Coeff of sine In 6^^^. 


Multiples of 


(Seconds o^ arc) 

1 

» 


F 

D 

r 

22639.5 

1 

0 

0 

0 

0 

-4586.426 

1 

0 

0 

-2. 

0 

2369.902 

0 

0 

0 

2 

0 

769.016 

2 

0 

0 

0 

0 

-668.111 

0 

1 

0 

.0 

0 

-411.608 " 

0 

0 

2 

0 

0 

4l 1.656 

2 

0 

0 

-2 

0 

-205.962 

1 

1 

0 

-2 

0 

-125.154 

0 

0 

0 

1 

0 

191.953 

1 

0 

0 

2 

0 

-165.145 

0 

1. 

0 

-2 

0 

147.693 

1 

1 

0 

0 

0 

-109.667 

1 

1 

0 

0 

0 

Coeff In sine in Latitude, 
(Seconds of arc) 

13461.48 

0 

0 

1 

0 

0 

1010.180 

1 

0 

1 

0 

0 

-999.695 

-1 

0 

1 

0 

0 

-623.658 

0 

0 

1 

-2 

0 

117.262 

0 

0 

1 

2 

0 

199.485 

-1 

0 

1 

2 

0 

-166.577 

1 

0 

1 

-2 

0 

61.913 

2 

0 

1 

0 

0 
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Multiples Of 
i t" F Dr 


f, X lO® of cosine In oo» 
and of sine In §in 


99972 

1674 

32 

1 

2 

-1675 

-4 

-2 

4 

-4 


0 10 0 1 

0 2 0 0 1 

0 3 0 0 1 

0 4 0 0 1 

0 10 11 

0 0 0 0 1 

0-1 0 01 

0 1 0-11 
0 0 1-10 
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The derived ecliptic lunar position is converted to Inertial 

coordinates 


V " aosJijBoa^^ 

n" ■ •inX.ooa 4_o<?ec - siM^ainc 
m m ’m 

V- ■ einijsoet + einXjsoai-einc 
^m m m 


where c is the obliquity to the ecliptic and 

« r + D + r + 
m m ♦ 

The conversion to earth-fixed coordinates is accomplished by a matrix 
transformation 



Inertial 


( 28 ) 


where 


aosQ^ einQ^ 0 

9 9 

M ® I -ei/iGg cos0g 0 


and 0 is the Greenwich hour angle (Table ii). 

A 

The unit vector in the direction of the sun is derived from 
Newcomb's theory in the same manner. The ecliptic elements, 1, F, D, r , 
are the same as for the lunar development. The solar coordinates 
{oosX »sin\.) are expressed as algebraic sums as listed in Table IV 

9 9 

The conversion to equatorial coordinates is 


1/ 
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X" • oosX^ 
s 

p" * sinX^ooae 
s 

v" = sinX sine 
«# 

and a final transformation to earth-fixed coordinates is applied, using 
Equation (28). 

A comparison of the radial body tide displacement computed by this 
analytic procedure and the ephemeris calculation from 6E0DYN is displayed in 
Figure 29. 
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4.0 NEW TECHNOLOGY 

The effort under this contract consisted of the development and pro- 
gramming of techniques to numerically calculate earth response to global 
semidiurnal and diurnal ocean tide models. Global Vertical crustal defor- 
mations have been evaluated for Mg, Sg, Ng, Kg, K^, 0^ , and ocean tide 
loading^ while horizontal deformations have been evaluated for the Mg tidal 
load. Tidal gravity calculations were performed for Mg and 0^ tidal loads, 
and strain tensor elements were evaluated for Mg loads. The Mg solution used 
for the ocean tide Included the effects of self-gravitation and crustal 
loading. 

Frequent reviews and a final survey for new technology were performed. 

It is believed that the mathematical and programming techniques and 
algorithms developed do not represent "reportable items," or patentable 
items, within the meaning of the New Technology Clause. Our reviews and 
final survey found no other items which could be considered reportable 
items under the New Technology Clause. 
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FIGURE 3: N 2 Tide Model neglecting ocean loading and self-gravitation 
effects Integrated on a 3*x3* global grid, 
g: Dashed lines — - cotidal lines (phase) In hours 

A: Solid lines corange lines (anplltude) in centiwters 
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Figure 8 

M 2 Vertical Displacement 

Dashed Lines cotidal lines (phase) In hours 
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Figure 12 
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* A(^,X) cos [o(t-to) - A + e] 
o =• .00007292 




>1 Vertical Displacement 

Dashed Lines cotidal lines (phase) in hours 

Solid Lines corange lines (amplitude) In aiUlneters 











Figure 14 

Vertical Displacement 

Dashed Lines cotidal lines (phase) in hours 

Solid Lines corange lines (amplitude) in millimeters 

U(4>,X;t) ■ A(^,X) cos [<Kt”To) - A + e] 
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Figure 16 

M 2 Horizontal Displacement in"! Direction 
Dashed Lines cotidal lines (phase) in hours 

Solid Lines corange lines (amplitude) in tenths of Billlaeters 

V, (i}>,X;t) = A(4i,X) cos [a(T-xo) - A + e] 
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Figure 22 

Horizontal Surface Strain Component over the Ocean 
a: Dashed Lines — cotidal lines (phase) in hours 
A: Solid Lines corange lines (amplitude) (xlO 








Figure 24 

Horizontal Surface Strain Component over the Ocean 














Figure 26 

Horizontal Surface Strain Component over the Ocean 
V: Dashed lines — cotidal lines (phase) in hours 












Global Horizontal Surface Strain Amplitude 
A: Solid Lines corange lines (amplitude)(xlO 
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